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Abstract. Using steady, axisymmetric, ideal magnetohydrodynamics (MHD) we 
analyze relativistic outflows by means of examining the momentum equation along 

' the flow and in the transfield direction. We argue that the asymptotic Lorentz factor 

is 7oo ~ /i — o"m, and the asymptotic value of the Poynting-to-matter energy flux 
ratio - the so-called a function - is given by <7oc/(l + ffoa) ~ itm/Mj where <tm is the 
(**] ' Michel's magnetization parameter and \ic 2 the total energy-to-mass flux ratio. We 

. discuss how these values depend on the conditions near the origin of the flow. By 

employing self-similar solutions we verify the above result, and show that a Poynting- 
00 \ dominated flow near the source reaches equipartition between Poynting and matter 

energy fluxes, or even becomes matter-dominated, depending on the value of ctm/m- 

t-H ■ 
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1. Introduction 

o ■ 

The main driving mechanism for relativistic outflows in AGNs, GRBs, 
and pulsar winds, is likely related to magnetic fields. These fields are 
^ \ able to tap the rotational energy of the disk, and accelerate matter 

5_i ■ ejecta not only magnetocentrifugally, but also due to the magnetic 

c/2 , pressure. Using ideal magnetohydrodynamics (MHD) we examine ini- 

tially Poynting-dominated outflows, trying to answer the following ba- 
sic question: Which part of the total ejected energy flux is transfered 
to the matter kinetic energy flux asymptotically, and how this value 
jjj ■ depends on the conditions near the origin of the flow? 

The system of equations of special relativistic, steady, cold, ideal 
MHD, consist of the Maxwell equations = V'B = Vx£ = Vx 
B - AirJ/c = V • E - iirJ°/c, the Ohm's law E = B x V/c, the 
continuity V • (po^yV) = 0, and momentum — 7Po (V • V) ('yV) — VP + 
(J°E + J x B) /c = equations. Here V is the velocity of the outflow, 
7 the associated Lorentz factor, (E , B) the electromagnetic field as 
measured in the central object's frame, J°/c , J the charge and current 
densities, and po the gas rest-mass density in the comoving frame. 

Assuming axisymmetry [d/dqb = 0, in cylindrical (z,uj,(fi) coordi- 
nates with z along the rotation axis] , four conserved quantities along the 
flow exist. If A = (l/27r) JJ B p ■ dS is the poloidal magnetic flux func- 
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tion, they are (e.g., Vlahakis & Konigl, 2003a) 1 : the mass-to-magnetic 
flux ratio ^a(A) = ^ir'ypoVp/Bp, the field angular velocity Q(A) = 
(Vtp/zu) — (V p /zu)(B ( f,/B p ), the specific angular momentum L(A) = 
^wV^ — zuB^/^a, and the energy-to-mass flux ratio 



n(A)c 2 = 7c 2 - mSlB^/^A ■ 



The right-hand side of eq. (1) consists of the matter energy-to-mass 
flux ratio 7c 2 , and the Poynting-to-mass flux ratio (fi — 7)c 2 . The sum 
of these two parts is a constant of motion, while their ratio is defined 
as the Poynting-to-matter energy flux ratio a = (fj, — 7) /j. 

All the physical quantities can be written as functions of (A, a): 
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where x = zu£l/c is the cylindrical distance in units of the light cylin- 
der's lever arm, and xa = (L£l/ pc 2 ) 1 / 2 its value at the Alfven point. 

The functions A(zu,z), a(zu,z) obey the two remaining equations 
of the system: The Bernoulli equation 2 which is a quartic for a 



l + a 



1 + 



q 2 w\va\ 



1 - 



+ <T- 



1 



3T 



+ 



,-y>2 

x A 



H- 



1+a 



,(5) 



and the transfield force-balance equation 
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Here the operator L = V 2 — ^zu ■ V is related to the curvature radius 
of the poloidal field lines TZ = \VA\ (LA — VA • V In \VAfzu\y 1 . 

1 The subscripts p/(j> denote poloidal/azimuthal components. 

2 This equation comes from the identity 7 2 = 1 + (7V p /c) 2 + ('yV ( p/c) 2 . 
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In the force-free limit a = oo it is \i = oo, \I> a = 0, L = oo, x\ = 1, 



while n/a, ji^ a, L^a, c(1 



are finite. In this case, only the first 



four terms of eq. (6) survive, resulting in the "pulsar equation" 

zuLA w ■ VA ^ac 6 d h^a vo\VA\ d£l 
\VA\ + \VA\ ~Wm\VA\ dA~tT + SI dA 



0. (7) 



2. The a function and its asymptotic value 

2.1. A GENERAL ANALYSIS 

An important combination of the field line constants is the "Michel's 
magnetization parameter" <tm(-<4) = AQ /^^c 3 . In terms of cm, and 
using eqs. (1) and (2), we may write the exact expression 3 



1+CT 




(8) 



The left-hand side represents the Poynting-to-total energy flux, and - 
using the first of eqs. (3) - can be rewritten as (/x — 7)//i- As long as 
the flow is Poynting-dominated (7 <C n), this ratio is close to unity, 
and a S> 1. This continues to be the case in the neighborhood of the 
classical fast magnetosound surface, where 7 w /i 1 / 3 (e.g., Camenzind, 
1986). As a result, the superfast regime of the flow is the only place 
where a transition from high (s> 1) to ~ 1 or <g 1 values of a is possible. 
In this regime, and for extremely relativistic flows, the term (—B^/E) 
is very close to unity. 4 Hence, eq. (8) gives a simple relation between 
the a function and the ratio B p w 2 /A: 

(9) 

Suppose that the value of the function B p w 2 /A near the classical fast 
surface is (B p zu 2 /A)f. Since a S> 1 at this point, eq. (9) implies that 
the constant of motion ctm/a* ~ 1/ '{B v w 2 / 'A) r. 

Denoting with 5l± the distance between two neighboring poloidal 
field lines A and A + 5 A, magnetic flux conservation implies B p w 2 = 
(w/5l±)5A. Thus, a decreasing B p w 2 - and hence, from eq. (9), a 

3 Eq. (8) remains the same with thermal effects included. 

4 The requirement that the Lorentz invariant B 2 — E 2 > 0, using E — xB p , gives 
Bl/E 2 > 1-1/x 2 . In addition, eqs. (3)-(4) give V^/c = x + (V p /c)(B^,/B p ), and the 
condition Vj, > implies -B^/E < c/V p . Thus, (1 - l/x 2 ) 1/2 < -B^/E < c/V p . 
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decreasing a function - corresponds to poloidal field lines expanding 
in a way such that their distance Slj_ increases faster than w. How 
fast the field lines expand is determined by the transfield force balance 
equation; thus, eq. (6) indirectly determines the flow acceleration. Since 
the available solid angle for expansion is finite, there is a minimum value 
of the B p w 2 /A function. The field lines asymptotically have a shape z ~ 
zq(A) +w/ tan-d(A), where $(A) is their opening angle. Differentiating 
the latter equation we get B p w 2 /A = (A-d'/siwd — Az' sinfl/w) , 
a decreasing function, reaching a minimum value sin-d/Ati' at w 3> 
z^sm 2 ®/!}' (Vlahakis, 2004). Since the factor sintf/Ai?' is ~ 1, the 
minimum value of the B p w 2 /A function is ~ 1, corresponding to 5 




(j?pCT 2 )oo a M 1 
(B p w 2 ) f ~ Li ~ (B p w 2 /A) f ■ 



(10) 



Equivalently, the asymptotic Lorentz factor is = n/(l + a^) ~ 
li — o"M) and the asymptotic Poynting-to-mass flux ratio is ~ omc 2 . 

Another interesting connection with the boundary conditions near 
the source can be found, by noting that, as long as jB^j E = xB p , 
B p w 2 /A rs 2|/|/.Af2, where |/| = (c/2)w\B ( p\ is the poloidal current. 
Thus, (B p w 2 /A) f ~ 2\I\f/A£l, and since |/| remains constant of mo- 
tion inside the force-free subfast regime, (B p w 2 /A) f « 2\I\i/A£l, and 
li/(Tm ~ 2\I\i/AQ. Hence, eq. (10) implies a direct connection of the 
efficiency and the Lorentz factor to the ejection characteristics 

croo An ( b v w 2 \ An ( An \ . . 



l + ^oo 21/L \ A / 21/1, ,uu ^ V 21 J„ 



2.2. The value of 0-qo in known solutions 

Solving the system of equations (5) and (6) is highly intractable. In the 
following we review the currently known methods to obtain solutions. 

2.2.1. Numerical methods 

Eq. (5) is a relatively simple algebraic equation for a. After substi- 
tuting its solution (in terms of A and its derivatives) in the transfield 

The only exception to this general result is to have asymptotically a/(l + a) <C 
(jm/h in some finite solid angle regions, combined with other regions with bunched 
field lines [in which B p zu 2 S> A and <r/(l + a) S> cm///]- Note also that the most 
general asymptotic field line shape slightly deviates from straight lines, resulting 
in a logarithmic acceleration reaching a values smaller than cxm/a* (Chiueh, Li, & 
Begelman, 1991; Okamoto, 2002; Vlahakis, 2004). However, this acceleration can 
happen in exponentially large distances and hence is physically irrelevant. 
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force-balance equation (6), we get a second order partial differential 
equation for the magnetic flux function A. Its solution determines the 
field-streamline shape on the poloidal plane. Due to the fact that this 
equation is of mixed type, i.e., changes from elliptic to hyperbolic, it is 
beyond the capability of existing numerical codes to solve this highly 
nonlinear problem, and no solution has been obtained so far. 

An alternative numerical approach is to solve the time-dependent 
problem (hyperbolic in time) and expect to reach a steady-state. How- 
ever, all existing codes fail to simulate relativistic magnetohydrody- 
namic flows for more than a few rotational periods. On top of that, it 
is not always clear how the issue of the boundary conditions is handled. 

A promising combination of the two above methods is followed by 
Bogovalov (2001), who solves the inner problem using time-dependent 
evolution (avoiding the elliptic to hyperbolic transitions) , and the outer 
problem using steady-state equations. The code is not yet capable of 
solving the problem at large distances, though. 

2.2.2. The force-free assumption 

In the force-free limit a = oo, the two equations (5) and (6) decouple. 
Thus, one may solve the somewhat simpler, elliptic equation (7) (e.g., 
Contopoulos et al., 1999), and then solve eq. (5) for a. 

The force-free solutions have several problems, the most important 
of which is the following: Since the back reaction of the matter to the 
field is neglected, the drift velocity soon after the light cylinder becomes 
larger than the light speed. Equivalently, there is no frame of reference 
where the electric field vanish, and eq. (5) has no real solutions for a. 

In any case, the force-free assumption brakes down in the superfast 
(7 > /i 1 / 3 ^ a < [i 2 / 3 — 1) regime, where the flow becomes hyperbolic 
and the back reaction of the matter to the field cannot be neglected. 
Since for Poynting-dominated flows the value of the a function at the 
classical fast surface is fi 2 ^ 3 -1 > 1, the force- free method cannot be 
used for examining the efficiency of the magnetic acceleration. 

2.2.3. The prescribed field line shape assumption 

If one assumes a known magnetic flux distribution, i.e., a known func- 
tion B p w 2 1 A = zu\VA\/A, then it is trivial to solve eq. (5) for a and 
find the flow speed 6 (e.g., Takahashi et al., 1990; Fendt & Greiner, 
2001). Thus, when we use this method, practically we implicitly give 
the function a\ However, these solutions do not satisfy the transfield 
force-balance equation (6); thus, they are not fully self-consistent. 



In the superfast regime, this equation reduces to the much simpler eq. [9]. 
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2.2.4. The monopole approximation 

This is a subcase of the prescribed field line case, based on the as- 
sumption that the poloidal magnetic field is quasi-monopolar. This 
assumption, with the help of eq. (9), is equivalent to the assumption 
that the a function is constant! In fact, a tiny acceleration is possible 
in the subfast regime, leading to ~ /U 1 / 3 and ~ /i 2 / 3 ^> 1 
(Michel, 1969). This solution gave the erroneous impression to the 
community that relativistic MHD is in general unable to give high 
acceleration efficiencies. However, the solution corresponds to a special 
case of boundary conditions, and most importantly, it does not satisfy 
the transfield force-balance equation. 

2.2.5. Other approximate solutions 

Various tries to solve a simplified version of eqs. (5) and (6), by neglect- 
ing some terms [e.g., Tomimatsu & Takahashi (2003) who neglected 
the curvature radius term in eq. (6)] involve a risk, since the system of 
equations is highly nonlinear and it could be crucial to keep terms that 
seem at first negligible, especially second order terms. 

Also perturbations around a monopolar solution [e.g., Beskin et al. 
(1998) who perturbed a solution with 7 = 00 and = 0] are (using 
eq. [9]) equivalent to perturbations of a around a constant high value. 

Concluding, in order to solve the efficiency problem, one has 
to solve simultaneously eqs. (5) and (6). 

2.2.6. The r self-similar special relativistic model 

The only known exact solution of eqs. (5) and (6) is the r self-similar 
special relativistic model, found independently by Li, Chiueh, & Begel- 
man (1992) and Contopoulos (1994) in the cold limit, and further 
generalized by Vlahakis & Konigl (2003a) including thermal effects. 
It corresponds to boundary conditions in a conical surface (6 = 9i 
in spherical coordinates [r, 6, (j)}) of the form B r = C\r F ~ 2 , = 
-C 2 r F ~\ V r = C 3 , V e = -C 4 , = C 5 , po = C 6 r 2 ^" 2 ), P = C 7 r 2 ^" 2 ), 
with constant C\,...fi-j. The parameter of the model F controls the 
initial current distribution [—wBj, = C2 sin 0i r F ~ l is an increasing or 
decreasing function of r for F > or < 1; see Vlahakis & Konigl (2003a) 
for details]. Despite the assumed form of the boundary conditions, the 
assumption that gravity is negligible, and the absence of intrinsic scale, 
r self-similar remain the only self-consistent relativistic MHD solutions. 

Figure 1 shows the function B p m 2 /A for various r self-similar so- 
lutions with application to GRB and AGN outflows. With the help of 
the figure and the known /x, h/gm values, we are able to fill the table 
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Figure 1. Function B p vo 2 /A versus 7 for exact r self-similar solutions presented in 
Vlahakis & Konigl (2003a, 2003b, 2004). The positions of the Alfven and classical 
fast magnetosound points are marked with "*" and "x", respectively. 



Table I. The efficiency in r self-similar models 



Solution 


/" 




/ B p zu 2 \ 




(B p 1jJ 2 )oo 




cru 


\ A )f 


\ / OO 




1 + <Joo 


Ia a 


10116.1 


2.02 


1.99 


1.49 


0.74 


0.73 


Ib a 


9997.4 


2.00 


1.99 


1.38 


0.69 


0.69 


Id a 


2150.0 


2.15 


2.06 


1.22 


0.59 


0.57 


IIa b 


778.9 


3.31 


3.22 


1.16 


0.36 


0.35 


IIb b 


1156.6 


3.06 


3.02 


1.13 


0.37 


0.37 


IIP 


75.0 


2.06 


1.93 


1.04 


0.54 


0.51 



a Trans- Alfvenic GRB outflow solutions, from Vlahakis & Konigl (2003a). 
6 Super- Alfvenic GRB outflow solutions, from Vlahakis & Konigl (2003b). 
c AGN jet solution, from Vlahakis & Konigl (2004). 



I. 7 This table verifies the results of § 2, and in particular: 1) Since the 
values of the 3rd and 4th columns are approximately equal, cjm/ 'fJ> ~ 
1 / {B p w 2 1 A) f indeed holds. 2) The value {B p w 2 / A)^ is indeed of order 
unity. 3) From the two last columns we verify that the asymptotic value 
of a satisfies eq. (9), and (Too /{I + ^oo) ~ (B p -cj 2 / A)J 1 (B p vd 2 /A)^ 
(<tm/ ^){BpW 2 /A)oo. 4) The value <Too/(l + cTqo) is roughly equal to 
<7m/m- 5) The value 700 is roughly equal to fj, — ctm = A*(l — cm/a 1 )- 
A common characteristic of all trans- Alfvenic r self-similar solutions is 
that /u/cjm ~ 2, and thus 700 ~ fj,/2. 

7 Besides the values of B p zu' 2 /A, fig. 1 shows 7^; we also find a x = (fl/foo) — 1. 
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2.2.7. z self-similar asymptotic solutions 

Vlahakis (2004) derived z self-similar asymptotic solutions of the sys- 
tem of eqs. (5) and (6) that verify the presented analysis. For con- 
ditions near the classical fast magnetosound surface corresponding to 
(B p w 2 ) /»v4<^>cr M </U<^> \I\i > AQ./2 the is found to be < 1. 

3. Conclusion 

• In order to solve for the acceleration it is absolutely necessary to 
solve for the poloidal field line shape as well. The Bernoulli and trans- 
field force-balance equations are interrelated and we cannot solve them 
separately, especially in the superfast regime. 

• Models that assume quasi-monopolar magnetic field (the Michel's so- 
lution included), equivalently assume that the magnetic acceleration is 
inefficient. Moreover, no fully consistent solution of this kind is known. 

• o"oo/(l + o-qo) ss (o"m/ 'fJ-)(B p w 2 1 '^4)oo ~ cm/m is an analytic expres- 
sion for the asymptotic Poynting-to-total energy flux ratio, confirmed 
by self-similar solutions. Possible acceleration to smaller can only 
happen in exponentially large (and thus physically irrelevant) distances. 

• Coo/(l + 0"oo) ps {AQ,/2\I\i)(B p w 2 / A)^ ~ AQ/2\I\i is a general result 
connecting the asymptotic value <Too (and hence = /Vfl + ^oo]) with 
the conditions near the origin of an initially Poynting-dominated flow. 

• The asymptotic Lorentz factor is 700 ~ [i — ayi ~ n(l — A£l/2\I\i), and 
the asymptotic Poynting-to-mass flux ratio is ~ ctmc 2 ~ fic 2 AQ/2\I\i. 
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